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Abstract 

Using data from the FOCUS experiment (FNAL-E831), we study the decay of 
A+ baryons into final states containing a A hyperon. The branching fractions of 
A+ into Air + , Ait + it + it~ and AK°K + relative to that into pK~n + are measured 
to be 0.217 ± 0.013 ± 0.020, 0.508 ± 0.024 ± 0.024 and 0.142 ± 0.018 ± 0.022, re- 
spectively. We also report new measurements of ^ A c +~ > A ^ = 1-09 ± 0-H ± 0.19, 

ffir?^+ +,r ? = 0.26±0.06±0.09 and r(Ae+ ^l 169 ?^ )Jf+) = 0.33±0.10±0.04. 
Further, an analysis of the subresonant structure for the A+ — > kir + ir + ir~ decay 
mode is presented. 



1 Introduction 



During the past several years there has been significant progress in the exper- 
imental study of hadronic decays of charmed baryons. However the precision 
on branching fraction measurements is only about 40% for many Cabibbo- 
favored modes and even worse for Cabibbo-suppressed decays [1]. As a result, 
we are not yet able to distinguish between the decay rate predictions made by 
different theoretical models, e.g., the quark model approach to non-leptonic 
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charm decays and the Heavy Quark Effective Theory (HQET) [2,3,4]. In this 
paper we present a study of A+ baryons produced by the FOCUS experiment. 
We present improved measurements of the branching fractions of the Cabibbo- 
favored decays A+ -> A7T+,A+ -> A7r + 7r+7r~ and A+ -> AK°K + . From the 
measurement of the first two modes, we are also able to extract the relative 
branching ratios of the two decays A+ — > £°7r + and A+ — > £°7r + 7r + 7r~. We re- 
port a new measurement of the subresonant mode A+ — > 2(1690)° K + . Finally 
we present the first study of the subresonant structure of the A+ — > A7r + 7r + 7r~ 
decay mode. 



2 Event Reconstruction 

This analysis uses data collected by the FOCUS experiment during the 1996- 
1997 fixed-target run at Fermilab. 

FOCUS is a photo-production experiment equipped with very precise vertex- 
ing and particle identification detectors. The vertexing system is composed of 
a silicon microstrip detector (TS) embedded in the BeO target segments [5], 
and a second system of twelve microstrip planes (SSD) downstream of the tar- 
get. Downstream of the SSD, five stations of multiwire proportional chambers 
and two large aperture dipole magnets complete the charged particle tracking 
and momentum measurement system. Three multicell threshold Cerenkov de- 
tectors are used to identify electrons, pions, kaons, and protons. The FOCUS 
apparatus also contains one hadronic and two electromagnetic calorimeters as 
well as two muon detectors. 

All decay modes reported have a A hyperon 1 in the final state. A detailed 
description of A and Kg reconstruction techniques in FOCUS is reported in 
Reference [6]. 

Candidates are reconstructed by first forming a vertex with tracks consistent 
with a specific A+ decay hypothesis. A cut on the confidence level that these 
tracks form a good vertex is applied. Production vertex candidates are found 
using a candidate driven vertexing algorithm which uses the A+ candidate 
momentum to define the line of the flight of the charm particle [7]. This 
seed track is intersected with other tracks in the event to form a production 
vertex. The confidence level for the production vertex must be greater than 
1%. Most of the background is rejected by applying a separation cut between 
the production and decay vertices: we require the significance of separation 
between the two vertices, L/a L , to be greater than some number, depending 

1 Throughout this paper the charged conjugate state is implied unless explicitly 
stated. 
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on the decay mode. 



All charged microstrip track segments from the charm decay must be linked 
to a single multi-wire proportional chamber track segment, be of good quality 
and be inconsistent with zero degree tracks from photon conversion. The like- 
lihood for each charged particle to be proton, kaon, pion or electron based on 
Cerenkov particle identification is used to make additional requirements [8]. 
For pion candidates, we require a loose cut that no alternative hypothesis is 
favored over the pion hypothesis by more than 6 units of log-likelihood. In 
addition, for each kaon candidate we require the negative log-likelihood kaon 
hypothesis, Wk = —2 In (kaon likelihood), to be favored over the correspond- 
ing pion hypothesis W w by W n — Wk > 3. 

The reconstructed mass of the A candidates must be between 1.1 and 1.125 
GeV/c 2 ; no cut is applied on the normalized mass [M(A) — M(A)pe>g\/< j m(A), 
because it is not centered around zero, probably due to the higher background 
under the signal region. We moreover require the higher momentum track used 
to reconstruct the A candidates to be compatible with the proton hypothesis, 
applying the cut W n — W p > 4. The reconstructed mass of the Kg must be 
within three standard deviations of the nominal Kg mass. 

We require A+ candidates to have a minimum momentum of 45 GeV/c and 
to have a lifetime less than five times the nominal value [1]. Finally, in order 
to reduce backgrounds, we require the production vertex to be located inside 
the target material. 



3 The normalization mode 



The A+ — > pK~ii + channel is our highest statistics A+ decay mode and it is 
used as the normalization mode for branching ratio measurements to mini- 
mize the overall statistical uncertainty. Moreover, all previous measurements 
in the literature [1] use this decay as a normalization mode, thus making any 
comparison straightforward. 

In order to minimize systematic biases, the normalization mode is selected us- 
ing the same cuts and the same fitting technique as the specific decay whenever 
possible. In addition, for each proton candidate we apply the cuts W n — W p > 4 
and Wk — W p > 1. The pK~ir + invariant mass distribution for an L/ctl > 4 
cut is shown in Fig. 1 (b). The resultant yield is 16447 ± 193 events. 
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Fig. 1. Invariant mass distributions for (a) A+ — > Att + and (b) A+ — > pAT 7r + . The 
fits are described in the text. 

4 The A+ — > An + decay mode 



We measure the branching ratio of A+ — > A7r + relative to A+ — > pK~7i + . In 
Fig. 1 (a) the An + invariant mass distribution for an L/o"l > 4 cut is presented. 
The confidence level for the decay vertex must be greater than 1%. We also 
apply a | cos#| < 0.6 cut, where 9 is the angle between the A momentum in 
the A+ rest frame and the A+ laboratory momentum. 

We note a broad structure around 2.2 GeV/c 2 coming from the decay mode 
A+ — > S°(A7)7r + where the photon from the S° decay is not reconstructed. 
The shape for this reflection has been obtained from a Monte Carlo simulation 
of this decay mode. The fit is performed using two Gaussians with the same 
mean for the signal, the reflection from the E°7r + mode, and a second order 
Chebychev polynomial for the background. The ratio of yields and the resolu- 
tions of the two Gaussians are fixed to the Monte Carlo values. The resultant 
yield is 750 ± 44 events. Correcting for the relative efficiencies estimated by 
our Monte Carlo simulation, we determine the branching ratio to be 



r(A+ - a 



7T 



r(A+ -> pK-ir- 



0.217 ±0.013 (stat. 



The number of fitted A+ — > S 7r + reflection events is 919 ± 92. Correcting for 
the relative efficiencies, we extract the relative branching ratio: 



T(A+ -> S°7r+) 

r(A+ - Att+) 



1.09 ±0.11 (stat.) 



(2) 
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Fig. 2. Invariant mass distributions for (a) A+ — > Att+7t+tt and (b) 
A+ — » Att+7t+7t~ for the subresonant analysis. The fits are described in the text. 

5 The A+ — > Att+tt+tt^ decay mode 



We measure the branching ratio of A+ — > Att+tt+tt^ relative to A+ — > pK~7i + . 
In Fig. 2 (a) the Att+tt+tt^ invariant mass distribution for an L/cr L > 5 cut 
is presented. The confidence level for the decay vertex must be greater than 
5%. We also apply a cos 9 > —0.9 cut, where 6 is the angle between the A 
momentum in the A+ rest frame and the A+ laboratory momentum. 

We also note in this decay mode a broad structure around 2.2 GeV/c 2 com- 
ing from the decay mode A+ — > £°(A7)tt+tt+tt- where the photon from the 
S° decay has not been reconstructed. This has been accounted for as in the 
A+ — > A?t + decay. The components of the fitting function are the same as 
in the A+ — > Att + case. The resultant A+ — > Att+tt+tt^ yield is 1356 ± 60 
events. Correcting for the relative efficiencies estimated by our Monte Carlo 
simulation, we determine the branching ratio to be 

T(A+ -> Att+tt+tt-) , 

r(k^ P K-* + ) = °- 5 ° 8 ± (stat °- (3) 



The number of fitted A+ — > X°tt+tt+tt reflection events is 480 ± 110. Cor- 
recting for the relative efficiencies, we extract the relative branching ratio: 

l^^KK 71 ! = 0-26 ± 0.06 (stat.). (4) 
T(A+ -> Att+tt+tt-) v ; v ; 



We have studied the subresonant structure in the decay mode A+ — > Att+tt+tt^. 
Considering our limited statistics, which would make a coherent analysis dif- 
ficult, we use an incoherent binned fit method [9] developed by the E687 
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Fig. 3. (a) Ait , (b) Att + and (c) ir + ir invariant mass distributions in the decay 
mode A+ Att+tt+tt". The yields are, respectively, 143 ± 27, 149 ± 28 and 317 ± 
68. 

Collaboration, which assumes the final state is an incoherent superposition of 
subresonant decay modes. 

For the resonant substructure analysis of A+ — > Att + 7t + 7t~ we enhance the 
signal to noise ratio applying an L/c^ > 8 cut and requiring 1.11 < M(A) < 
1.119 GeV/c 2 . In Fig. 2 (b) the A7r + 7r + 7r _ invariant mass distribution for 
events which satisfy these cuts is presented. The resultant yield is 594 ±31 
events. 

A study of the two-body invariant mass distributions was done to better iden- 
tify which resonances may contribute to the A7r + 7r + 7r~ decay channel. In Fig. 
3 the two body Ayr", Att + and tt + tt~ invariant mass distributions provide ev- 
idence for the E^SS) 1 * 1 and p(770)° resonances. For this study we require 
the A7r + 7r + 7r~ invariant mass to be within 2a (18 MeV/c 2 ) of the A+ nominal 
mass and we perform a sideband subtraction to reduce the background. The 
fits are performed using Breit-Wigners for the signal shape, with the mean 
and width fixed to the Monte Carlo values, and Chebychev polynomials for 
the backgrounds. 

For subresonant modes in the resonant analysis we therefore consider the 
channels £(1385)-7r + 7r + , £(1385)+7r+7r-, Ap(770)°7r+ and £(1385)+p(770)°, 
plus a nonresonant channel (A7r + 7r + 7r _ )jvi?- All states not explicitly considered 
are assumed to be included in the nonresonant channel. 

We determine the acceptance corrected yield into each subresonant mode using 
a weighting technique whereby each event is weighted by its kinematic values in 
the three submasses (Ati~), (Att + ) and (7r + 7r~). We construct eight population 
bins depending on whether each of the three submasses falls into the expected 
resonance peak (within the nominal width). From a Monte Carlo simulation 
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Table 1 

Fractions relative to the inclusive mode for the subresonant structure of the A+ 
A7r + 7r + 7r~ decay mode. 



Subresonant Mode 


Fraction of — > A7r^7r^7r 


(A7r + 7r + 7r~)jvii 
Air + p 


< 0.30 @90% CL 
0.21 ± 0.03 ± 0.02 
0.28 ± 0.10 ± 0.08 
0.40 ± 0.12 ± 0.12 
0.14 ± 0.09 ± 0.07 



of each subresonant mode a, we compute the bin population rii in the eight 
bins and we calculate a transport matrix T ia between the number of generated 
Monte Carlo events Y a and the bin populations: 

n i = ^ Tj a Y a . (5) 

a 



The elements of the T matrix can be summed to give the efficiency e a for each 
mode: 

£a=J2 T ia- (6) 



This Monte Carlo determined matrix is inverted to create a new weighting 
matrix which multiplies the bin populations to produce efficiency corrected 
yields. Each data event can then be weighted according to its values in the 
submass bins. Once the weighted distributions for each of the five modes 
have been generated, we determine the acceptance corrected yields by fitting 
the distributions with two Gaussians with the same mean and a second or- 
der Chebychev polynomial for the background. Using incoherent Monte Carlo 
mixtures of the five subresonant modes we verify that the method is able to 
correctly reproduce the generated mixtures of the different modes. 

The results for the A7r + 7r + 7r~ decay are summarized in Table 1. The five 
weighted histograms are shown in Fig. 4, where Fig. 4 (f) is the weighted dis- 
tribution for the sum of all subresonant modes. The systematic uncertainty 
for the subresonant fractions is estimated varying the width of the resonance 
peaks in the construction of the kinematic bins. The goodness of fit is evalu- 
ated by calculating a x 2 f° r the hypothesis of consistency between the model 
predictions and the observed data yields in each of the 8 submass bins. We 
obtain a x 2 of 7.86 (for 3 degrees of freedom) and a confidence level of about 
5%. 
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Fig. 4. A+ — > A7r + 7r + 7r weighted invariant mass distributions for (a) 
(Ait + it + it-) nr , (b) S(1385)-^+7r+, (c) S(1385)+7r+7r-, (d) A / 9(770)°7r+, (e) 
S(1385)+p(770)°, (f) inclusive sum of all five modes. 

6 The A+ A~K°K + decay mode 



We measure the branching ratio of A+ — > AK°K + relative to A+ — > pK~7i + . 
The K° are detected through K^s. Due to the limited phase space, the signal 
can be observed without the need for a L/ai or decay vertex confidence level 
cut. In Fig. 5 the AKgK + invariant mass distribution is presented. 

The fit is performed using two Gaussians with the same mean for the signal and 
a second order Chebychev polynomial for the background. The ratio of yields 
and the resolutions of the two Gaussians are fixed to the Monte Carlo values. 
The resultant yield is 251 ± 31 events. Correcting for the relative efficiencies 
estimated by our Monte Carlo simulation, we determine the branching ratio 
to be 

T(A+ -> AK°K+) . 

kk^vK^) ^ =0.142±0.018 (Stat.). (7) 
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Fig. 5. Invariant mass distribution for A+ — > AXgX + . The fit is described in the 
text. 

The Belle collaboration [10] has recently shown evidence of the resonant con- 
tribution A+ -> E(1Q90)°K + in the decay A+ -> Ai^iC+ with the 5(1690)° 
reconstructed in AiT°. In our analysis, the H(1690) i^ + events are selected 
using the same cuts used for the AKgK + mode; the AK$K + invariant mass 
is required to be within 2a (10 MeV/c 2 ) of the A+ nominal mass. A sideband 
subtraction is performed to reduce the combinatoric background under the A+ 
signal region. 

The AKg invariant mass distribution is shown in Fig. 6. The fit is performed 
using a Breit-Wigner function for the signal and a first order Chebychev poly- 
nomial for the background. The mean and the width of the Breit-Wigner are 
fixed to the Monte Carlo values. 2 The resultant yield is 84 ± 24 events. 

We measure the branching ratio relative to A+ — > AK°K + to be 

T(A+ -> S(1690)°K+) 



T(A+ -> AK°K-< 



x S(H(1690) U -> AK ) = 0.33 ± 0.10 (stat.). (8) 



7 Systematic studies 



The systematic effects are evaluated after investigation of different sources: 
uncertainties in the reconstruction efficiency and in the resonant substructure 
for multibody decays and the choice of fitting conditions. 

To determine the systematic error due to the reconstruction efficiency we 

2 The E(1690)° is generated in our Monte Carlo simulation with a mass of 1.688 
GeV/c 2 and a width of 10 MeV/c 2 . 
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Fig. 6. Invariant mass distribution for A.K S in the decay — ► KK S K + . The fit is 
described in the text. 



follow a procedure based on the S-factor method used by the Particle Data 
Group [1]. For each mode we split the data sample into independent sub- 
samples based on A+ momentum, data-taking period, particle-antiparticle, 
significance of separation between production and decay vertices and differ- 
ent A and Kg categories, based on the location and geometry of the neutral 
particle decay. These splits provide a check on the Monte Carlo simulation of 
charm production, of the vertex detector and of different variables employed 
in the event selection. We define the split sample variance as the difference 
between the scaled variance and the statistical variance if the former exceeds 
the latter. The method is described in detail in Reference [11]. 

Considering the large uncertainty on the measured subresonant fractions in 
the multibody decays, we also vary these fractions in the Monte Carlo simula- 
tion and we use the variance in the branching ratios as a contribution to the 
systematic error. 

We measure the systematic uncertainty due to fitting conditions using a fit 
variation technique, which includes variations in bin size, fitting range, back- 
ground and signal shapes (different order of the Chebychev polynomial, leaving 
the two Gaussian parameters free in the fit or using a single Gaussian for the 
signal) . 

We also include a systematic error contribution from the absolute tracking 
efficiency for the different multiplicities in the final states. In Table 2 we sum- 
marize the systematic uncertainty for each mode. Several measurements for 
the modes reported here are present in the literature [12,13,14,15,16,17,18]. 
In Table 3 we present the FOCUS results with a comparison to the PDG 
values [1]. 
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Table 2 

The systematic uncertainties from the Monte Carlo simulation, the fitting condition 
and the the total for each mode. 



Mode 


Simulation 


Subresonances 


Tracking 


Fit 


Total 


T~V A ~t" A + \ 

r(A+-^pif-7r+) 
r(A+^s°7r+) 
r(A+^A7r+) 

r(A+^A7T+7T + 7r-) 

r(A+^ P K-n+) 

r(A+^E°7r+7r+7r-) 
r(A+^A7r+7r+7r-) 
r(A+^AK°K+) 

r(A+^S(1690)°(AX°)i<r+) 


0.017 
0.19 
0.016 
0.08 
0.021 


— 
0.010 
0.001 


0.005 
0.004 


0.008 
0.04 
0.014 
0.03 
0.005 
0.04 


0.020 
0.19 
0.024 
0.09 
0.022 
0.04 


r(A+^AK°K+) 









8 Conclusions 



We have investigated and measured the branching ratios of several A+ Cabibbo- 
favored decay modes containing the A hyperon in the final state. These modes 
are A+ -> Att+, A+ -> Att+tt+tt- and A+ -> AK°K + . From the fit to the first 
two modes, we are also able to extract the relative branching ratios of the 
two decays A+ — > S°7r + and A+ — > £°7r + 7r + 7r~. These measurements are an 
improvement over previous results for the same decay modes. We report a new 
measurement of the subresonant mode A+ — > S(1690)° K + consistent with the 
recent Belle result. We have also performed an analysis of the subresonant 
structure of the decay A+ — > A7r + 7r + 7r~. We observe a small nonresonant com- 
ponent and the presence of vector resonances in the dominant modes, as it 
has been observed in most charm meson decays. 
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Table 3. FOCUS results compared to previous measurements. No direct measurement exists for the relative branching ratios 
r(A+ -> S°7r+)/r(A+ -> Avr+) and T(A+ -> S°7r + 7r+7r")/r(A+ -> Att+vt+tt"). The relative efficiency includes the branching frac- 
tions into the observed final state particles. 



A+ Decay 


Signal 


A+ Reference 


Reference 


Relative 


FOCUS 


PDG [1] 


Mode 


Yield 


Mode 


Yield 


Efficiency 






Att+ 


750 ± 44 


pK~n + 


16447 ± 193 


0.209 ± 0.001 


0.217 ± 0.013 ± 0.020 


0.180 ± 0.032 




919 ± 92 


Att+ 


750 ± 44 


1.119 ± 0.001 


1.09 ± 0.11 ± 0.19 


1.11 ± 0.49 


A7r + 7r + 7r _ 


1356 ± 60 


pK~7T + 


12898 ± 147 


0.207 ± 0.001 


0.508 ± 0.024 ± 0.024 


0.66 ± 0.11 


S°7r + 7T + 7r- 


480 ± 110 


A/T + 7r + 7r~ 


1356 ± 60 


1.375 ± 0.001 


0.26 ± 0.06 ± 0.09 


0.33 ± 0.16 


AK°K+ 


251 ± 31 


pK~7Y + 


10952 ± 132 


0.161 ± 0.001 


0.142 ± 0.018 ± 0.022 


0.12 ± 0.02 ± 0.02 


E{1690)°(AK°)K + 


84 ± 24 


AK°K + 


251 ± 31 


1.00 


0.33 ± 0.10 ± 0.04 


0.26 ± 0.08 ± 0.03 



